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Abstract 

We deal with the periodic boundary value problem for a second-order nonlinear 
ODE which includes the case of the Nagumo type equation v xx — gv + n{x)F(v) = 0, 
previously considered by Grindrod and Sleeman and by Chen and Bell in the study 
of nerve fiber models. In some recent works we discussed the case of nonexistence 
of nontrivial solutions as well as the case in which many positive periodic solutions 
may arise, the different situations depending by threshold parameters related to 
the weight function n(x). Here we show that for a step function n{x) (or for small 
perturbations of it) it is possible to obtain infinitely many periodic solutions and 
chaotic dynamics, due to the presence of a topological horseshoe (according to 
Kennedy and Yorke). 
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1 Introduction and statement of the main results 



This paper deals with the study of a class of nonlinear second order ODEs 
which was proposed by Chen and Bell [3] as a model for a nerve fiber with 
spines. In their model, the authors consider a nerve fiber as an infinitely long 
cable. The spines, which are small protrusions on the membrane, are assumed 
to be periodically distributed. Their population is characterized by a density 
n = n{x) and by a lumped ohmic resistance 1/ p of the spine stem. The current 
flowing through the spine stem is given by p(u — v), where v (x, t) and u(x,t) 
represent two membrane potentials (v is the one corresponding to the dendrite, 
while u is that at the spine head). Chen and Bell, following the cable theory 
approach previously considered by Baer and Rinzel [1] and simplifying the 
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active membrane dynamics from the Ho dgkin- Huxley form to the FitzHugh- 
Nagumo's, derived the following system 



v t = v xx - gv + n(x)p(u - v) 
u t = p(u -v) + f(u), 

where g > is a constant such that 1/g represents the dendritic membrane 
resistance and / has the typical cubic form of the Nagumo equation [13] . The 
search of the steady state solutions u = u(x),v = v(x) of the above system, 
lead to the study of 

v X x ~ gv + n(x)f(u) = 
v = S(u;p) :=u- 

In [3], the authors discuss, separately, the cases in which the function S is 
invertible or not (for a given p). If S is invertible and therefore we can write 
u = R(v), for a suitable monotone increasing function R, we can confine 
ourselves to the study of the nonlinear second order scalar ODE 

Vxx - gv + n(x)F(v) = 0, (1) 

with F(v) := f(R(v)). We note that the possibility of inverting S and therefore 
considering equation (1) is not too restrictive and can be assumed when p is 
sufficiently large. In [3] the function n(x) is supposed to be piecewise constant 
and /3-periodic. More precisely, two positive constants n (small) and n\ (large) 
are taken such that 



n(x) = < 



ni , for < x < a, 
n , for a < x < f3. 



In [3, Lemma 4.1], assuming further a condition which implies that Jq 1 F(s) ds > 
0, it is stated that equation (1) has a positive periodic solution when a is suf- 
ficiently large (i.e., a is close to (5). 

In some recent papers [20,21] we extended the above recalled Chen and Bell 
result toward different directions. In particular, in [20], we proved (via a vari- 
ational approach) that the existence of at least two positive periodic solutions 
is guaranteed for an arbitrary weight function n(x), provided that Jq n{x) dx 
is sufficiently large. On the other hand, in [21], we obtained (via the Poincare- 
Birkhoff fixed point theorem) the existence of many positive periodic solutions 
under more restrictive hypotheses on n(x) (but general enough to include the 
piecewise constant case) and removing the requirement that Jq 1 F(s) ds > 0. 
In this work we come back to the original assumptions on n(x) of Chen and 
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Bell's paper with the aim to prove that their model presents solutions with a 
very complicated behavior. 



In order to introduce our results, we consider equation (1), where, throughout 
the paper, g > 0, n : R — > R + := [0, +oo) is a /5-periodic locally integrable 
function and F : K. — > K. is a C 1 function with a iV-shaped graph. More 
precisely, we assume that there exists a e ]0, 1[ such that 



and, moreover, 

F(s) < for < s < a and s > 1, F(s) > for s < and a < s < 1. 

A typical function F satisfying such conditions is given by the cubic nonlin- 
earity 



of the celebrated FitzHugh-Nagumo equation [8,9]. We look for solutions v(x) 
of (1) with 



and, in particular, we shall focus our attention to the search of periodic solu- 
tions as well as to the detection of solutions presenting some kind of chaotic- 
like behavior. 

Starting with the 70's a great deal of theoretical and numerical results have 
been concerned the investigation of nonlinear ordinary and partial differential 
equations and systems associated to the FitzHugh-Nagumo's. Variational or 
topological functional analytic methods have been successfully applied in the 
study of various boundary value problems associated to such equations. With 
this respect, we also recall some recent works by Sweers-Troy [18] and Dancer- 
Yan [4,5] where an abundance of solutions with a peculiar behavior for the 
Dirichlet or Neumann BVPs is established for related (although different from 
ours) Nagumo type systems. We stress, however, that some special features of 
equation (1), like the presence of a periodic weight, makes natural to address 
our investigation toward different aspects that, perhaps, have not been yet 
fully discussed in the literature. 

In the present paper we propose a dynamical system approach. To this end, 
we think at the one-dimensional re-variable (representing the longitudinal axial 
dimension of the idealized nerve fiber) as a temporal variable, in order to deal 
with a first order system in the plane of the form 



F(0) = F(a) = F(l) = 



F(s) = s(s-a)(l-s) 



0<v(x)<l, WxeR 



v 



y 



< 



(2) 




gv — n(t)F(v). 
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We denote by the Poincare's operator related to system (2), that is the map 
which associates to a point z G M 2 the value ip(z) := (((3,z) of the solution 
C(') — C(', z ) °f system (2) with £(0) = z. We notice that, without further 
assumptions on F, we don't have tp necessarily defined on the whole plane M 2 . 
On the other hand, by the smoothness of F, we know that ip is a (orientation- 
preserving) homeomorphism of its domain C M 2 onto its image i^CI 2 . 

Our main result is the following: 

Theorem 1.1 Assume that F(-) is strictly convex on [0, a], strictly concave 
on [a, 1] and satisfies 

[ F(s)ds>0. (3) 
Jo 

Define 

Let < a < (3 and suppose that n(-) is the stepwise function given by 



n(t) : = < 



rii , for < t < a, 

' 7 (5) 

n , /or a <t < (3, 



with < n < Wi . T/ien i/iere exzsfo c<j* > mj snc/i £/ia£ /or every n and n\ 
with 

< n < < to* < ni , 
t/iere exisi a compact invariant set 

A C]a, l[x]0,+oo) 

/or which %p\\ is semiconjugate via a continuous surjection h to a two-sided 
Bernoulli shift on two symbols. Moreover, for each periodic sequence of two 
symbols {ik)kei there is a point z G A which is periodic and such that h(z) = 
{ik)kez- 

As we shall see at the end of the proof of Theorem 1.1 (which is performed 
in Section 4) the semi conjugation of -?/>|a allows the following interpretation in 
terms of the solutions to (1): 

Consider a two-sided sequence of two symbols £ := {ik)k& with 

i k e{l,2}, WkeZ. 

Then, there is at least one solution = v^(x) of (1) which is defined for every 
i£l and satisfies 

< vz(x) < 1, ViGl 

as well as (v^(0), ^f^(O)) G A. Moreover, the symbol — 1 means that v^(x) 
has precisely two strict maximum points separated by one strict minimum point 
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along the time interval [{k — l)/3,(k — l)/3 + a], while, the symbol ik = 2 
means that v%(x) has precisely three strict maximum points separated by two 
strict minimum points along the time interval [{k — l)/3,(k — l)/3 + a]. In 
both the situations, v^(x) is convex in the interval [{k — 1)(5 + cx,k/3], with 
v'^((k — 1)(3 + a) < and v f ^(kj3) > 0. If the sequence (ik)kez is periodic, that 
is ik = ik+e for some £ > 1, then we can take v^(-) as a £{3-periodic solution 
as well. 

A minor modification in the argument of the proof of Theorem 1.1 yields to 
the following result. 

Theorem 1.2 Assume that F(-) is strictly convex on [0,a], strictly concave 
on [a, 1] and satisfies (3). Let be like in (4) and set 

9 

Suppose that n(-) is defined as in (5) with 

< n < m* < ml < n\ . 
Then there exist a* > and 5 > such that for every a and f3 with 

a*<a<P<a + 5, 
the same conclusion of Theorem 1.1 holds. 

We observe that it is possible to obtain extensions of Theorem 1.1 and Theo- 
rem 1.2 by producing symbolic dynamics on p objects with p > 2 (see Remark 
4.1). 



2 Notation and basic tools 

Throughout the paper we denote by R, R + and Rq the sets of real, real 
nonnegative and positive numbers. By 1 1 • 1 1 me mean a given norm (for instance, 
the euclidean one) in R 2 . 

In the sequel we consider some special subsets of the plane, called generalized 
rectangles. By such a name we mean any subset of R 2 which is homeomorphic 
to the unit square [0, l] 2 . It is possible to define an orientation for a generalized 
rectangle 1Z, by selecting two disjoint compact sub-arcs TZ7 and 1Z~ of its 
boundary, called the left and the right sides of 1Z. In a more formal way, 
using the Jordan-Shoenflies theorem, we can choose a homeomorphism h of 
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the plane onto itself such that 

h([0, l] 2 ) = 72, h({0} x [0, 1]) = 11J and h({l} x [0, 1]) = H~ . 
In this situation, we also set 

n~ := nj u n~ . 

By a path 7 we mean a continuous map 7 : [0, 1] — > R 2 . 

Let (.A, ^4~) and (B,B~) be oriented rectangles and suppose that ip : R 2 D 
— > R 2 is a map which is continuous on a set P C D^, n ^4. We say that 
(T>,ip) stretches (A, A~) to (B,B~) along the paths and write 

(V,^):(A,A-)^(B,B-), 

if, for every path 7 : [0, 1] — > A with 7(0) G .Af and 7(1) G „A~, there exists a 
subinterval [ti , t 2 ] C [0, 1] such that 

7 (0 eP, ^( 7 (t))eB, Vte[t u t 2 ] 

and -^(7(^1)), "^(7(^2)) belong to different components of £>~. 

It is easy to check that this property of stretching along the paths is preserved 
under the compositions of maps. Moreover, as proved in [15,16], if 

then there exists at least one fixed point for if) in T>. Such a fixed point property, 
when applied to different subsets of the domain and to the iterates of ip, may 
be exploited in order to find many different periodic points for the map ip. 

The key tool that we use for the proof of the existence of chaotic-like dynam- 
ics for the Poincare's map associated to equation (1) is the following lemma 
which combines the above mentioned fixed point theorem with the theory of 
topological horseshoes developed by Kennedy and Yorke in [11] (see also [10]). 

Lemma 2.1 Let (72-, 72") be an oriented rectangle and let ip : R 2 D — > R 2 
be a map. Assume there exist compact sets T>i C 72 fl (for i = 1,2) with 

V l r\V 2 = 0, 

such that 

(7^,V) : (72,72-) (72, TZ~) , Vi = 1,2. 
Then the following conclusions hold: 

• For any two— sided sequence of two symbols (ik)k& G S 2 := {1,2} Z , there 
exists a sequence (wk) keZ such that w k G V ik and ip(w k ) = w k +\ , for all 
keZ; 
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• // the sequence (ik)kez £ S 2 is m-periodic (m > 1), then there exists a 
corresponding sequence (wk)kez w ^h Wk G V ik which is m-periodic. 

Furthermore, as a consequence of the above properties, there exists a nonempty 
compact set 

ACV 1 UV 2 

which is invariant for ip and such that ?/>|a is semiconjugate to the two-sided 
Bernoulli shift on two symbols. The subset of A made by the periodic points of 
ip is dense in A and the counterimage (by the semiconjugacy) of any periodic 
sequence in S 2 contains a periodic point of ip. 

For the proof, see [15,16] and [17] for more recent details. Note that the semi- 
conjugation to the Bernoulli shift and the density of periodic points are typ- 
ical requirements for chaotic dynamics. The kind of chaotic-like dynamics 
described by Lemma 2.1 is based on the definition of chaos in the coin-tossing 
sense by Kirchgraber and Stoffer [12]. The same chaotic behavior is also ob- 
tained by Kennedy, Kogak e Yorke in [10, Proposition 5]. With respect to [10] 
and [12], our case takes into account also of the presence of periodic itineraries 
generated by periodic points. Concerning the applications to second order 
ODEs with periodic coefficients, where tp is the Poincare map associated to 
an equivalent first order differential system in the plane, the complex dynamics 
we obtain is in line with similar results appeared in the literature for different 
type of equations (see, for instance Capietto, Dambrosio and Papini [2]). Fi- 
nally, we refer to Mischaikow and Mrozek [14] and to Zgliczyhski and Gidea 
[22] for related topological approaches in higher dimension. 



3 Technical lemmas 



Let us consider the equation 



dx 

x" - gx + n(t)F(x) = (x' := — ), 
as well as the associated first order system in the phase-plane 



(6) 



x =y 

y' = gx — n(t)F(x), 



(7) 



where n(t) > is a /3-periodic piecewise continuous function that will be 
defined in a more precise manner in the sequel. 



The map F : 



is a C 1 function with a iV-shaped graph. In particular, 
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we assume there exists a G ]0, 1[ such that 



F(0) = F(a) = F(l) = 



and, moreover, 

F(s) < for < s < a and s > 1, F(s) > for s < and a < s < 1. 
A typical example for F is given by the Nagumo type cubic nonlinearity 



For the moment and in order to simplify the subsequent discussion, we assume 
(8). We point out, however, that the properties we are going to present below 
are still true for a broader class of functions. 

We are looking for solutions x(-) to (6) such that 



Accordingly, these solutions are the same also for any other equation where 
the function F(s) is the same like the given one in [0, 1] but possibly different 
elsewhere. For convenience we then suppose that 



so that all the solutions of equation (6) are globally defined. Clearly, at the 
conclusion of our argument, we need to check that the solutions we are inter- 
ested in satisfy the constraint in (9). To this purpose, we state the following 
lemma whose proof (which is valid also for a more general class of equations) 
is postponed in the Appendix. 

Lemma 3.1 Let n : [t , ^i] - ► R + be a Lebesgue integrable function and sup- 
pose that F : R — > M. is a locally Lipschitz function such that 



F (s) > for s < 0, F (s) < for s > 1 and F (s) = F(s), Vs G [0, 1]. 



F(s) = s(s - a)(l - s). 




< x(t) < 1, Vt e R. 



(9) 



F(s)\ < 1, Vs G R\ [0,1] 



Let x(-) be a solution of 



x" - gx + n(t)F (x) = 0, (g > 0) 



defined on [to, h] and such that 



< x(i ),a;(ii) < 1. 
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As a preliminary discussion, we start by a phase-plane analysis of the trajec- 
tories of (7) in the case when n(-) = constant = /i > 0. Accordingly, we study 
the autonomous system 



(E) 



x' = y 



y' = gx - /iF(x), 



and look at the effects of the parameter \i on the qualitative behavior of the 
orbits. 



First of all we observe that there is m* such that, for fi £ ]0, m* [ 

gs-nF{s)>0, Vs>0, 

while for each /z > there are two nontrivial equilibrium points P = (a M , 0) 
and Q = (c M , 0) to system (E) with 

0<a<a M <c M <l 

(see, for instance, Fig. 1 for an example which depicts this situation). More- 
over, a M — > a + and c M — > 1 _ as ^ — > +oo. The point P is a center and Q is a 
saddle. This is true for F(s) defined in (8) as well as for any function F which 
is strictly concave in [a, 1]. 

The actual value of mj can be computed as 

* _ 9 
m ° ~ F'(s*) ' 

where s* e ]a, 1[ solves the equation sF'(s) = F(s) yielding the maximal slope. 



System (E) is of conservative type with energy 

£(x, y) = £"(x, y) := \y 2 - g\x 2 + ^(x), (10) 

where x 

F{x) := [ F(s) ds. 
Jo 

We denote by L c the part of the energy line at level c lying in the strip 
< x < 1, that is 

L c ■■= {(x,y) : £(x,y) = c, < x < 1}. 
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Fig. 1. Different shapes for the graph of the function y = —gs + uF(s), for a 
small and a large value of [i. In this example, we have chosen g = 0.5 and 
F(x) = x(l — x)(x — a), with a = 0.4. The graph of the decreasing function corre- 
sponds to the case fx = no = 0.8, while the other graph (with two humps) is obtained 
for fi = n\ = 16. 

We present now some technical results which provide dynamical information 
about the trajectories of system (E) for different values of the parameter /i. 

Lemma 3.2 Suppose that /i e]0, mj[. Then, for every w = (po,0) with < 
Po < 1 the level line L c (c :— £{wq) ) passing through wq intersects any vertical 
line x = £ with po < £ < 1 exactly at two points ((,,—y^), (£, y^)- For < 
Po < 1, L c is an orbit path and the time needed to move from (£, — y^) to wq 
(which coincides to the time needed to move from w to (£, y^) ) along the orbit 
is given by 



o-(po,0 = I 



ds 



Jg{s*-ti)-2rtr(s)-f(p ))' 



Assume that 



(Ho) 
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Fig. 2. Level lines (at different energy) of system (E) with u = 10 (the closed lines) 
overlapped with those of the same system for u = 0.1. For this example, we have 
chosen g = 0.1 and F(x) = x(l — x)(x — a), with a = 0.4. 



Then, there is ml > m^ such that, for each /i > m* , the level line L is a 
homoclinic orbit of system (E) intersecting the abscissa at a point 0) with 



Moreover, for each c g]£(P),0[, the level line L c is a closed curve which 
corresponds to a periodic orbit of system (7) having minimal period 

2 ds 
= 7^ L ^2(F(x Q )-Hs))-l(xl-s 2 ) ' 

where xq and xi, satisfying 

< x < a M < xi < 6 M , 

are the two solutions of the equation £(x, 0) = c with x € ]0, 1[ . 

Proof. We start by considering system (E) for /i e]0, m^. Let L c be a level 
line passing through a point wo = (po,0) with < po < 1. By definition, we 



have 

The function Q is strictly increasing on [0, 1]. In fact G'(x) = 2(gx — fiF(x)) > 



y 2 = G(x) :=g(x 2 - Po 2 )-2ti(F(x)-F(p )), xe [0,1]. 



0, Vx G]0, 1]. Since G(po) — 0, it turns out that yG(x) is defined on [po, 1] 
and therefore L c intersects any vertical line x = £ (for po < £ < 1) exactly at 



the two points (£, — y^) and (£, with = y (?(£). From the relation 

which is satisfied for any time interval I in which y(t) = x'(t) > 0, we can 
compute the time cr(p ,£) at which a solution of system (E) (starting at the 
point w for t — 0) reaches the point (£, Actually, this is the standard 
time-mapping formula which reads in our case as 



G(x(t)) J* few' 

This proves the first part of our lemma. 

We consider now system (i?) for \i > uiq . In this case, we have two equilibrium 
points P = (a M , 0) and Q = (c M ,0), with a < a M < c M < 1. The level line L 
passing through the origin is defined by the relation 

y 2 = H(x) :=gx 2 - 2fiJ r (x), x e [0, 1]. 

The function 7i is strictly increasing and positive on ]0, a M ], strictly decreasing 
on [a M ,c M ] and strictly increasing on [c M , 1]. Hence, 7^(1) — g — 2pjF(l) > 
gc^ — 2pjF(c M ) = 7i(c M ). Since .F(l) = Jq F(s)ds > (by assumption (if ) ) 
we know that for /i sufficiently large, say fi > > ^ ne f unc ti° n 'H has exactly 
one zero in ]0, c M [ and such zero, denoted by 6 M , is contained in ]a M , c M [ . Thus, 
if we define as 

1- 2/ 1 F(a)da' 

we have that for every /j, > m*, the level line L is a close curve (the locus of 
a homoclinic trajectory of (E) ) intersecting the abscissa at the point (6 M ,0), 
with a M < 5 M < c M . A similar argument shows that, for every c e]£(P),0[, 
the level line L c is a closed curve surrounding the equilibrium point P and, as 
a matter of fact, it is the locus of a periodic orbit. 
In this case, the energy relation reads as 

y 2 = 2c + H(x), are [0,1]. 

The equation H(x) = —2c, with c e]£(P),0[ and x e]0, 1[ has exactly two 
solutions, xo and x\ , with 

< x < a M < xi < bf, . 
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We rewrite now the energy relation as 

y 2 = -H{x ) + H{x) = /i(2(^(x ) - F{x)) - ^(x 2 - x 2 )). 

Arguing like in the first part of this proof (the case of [i small) we can compute 
again the time-mapping related to the orbit L c . We find that the fundamental 
period t(x ) of such orbit is twice the time to move along L c from (x , 0) to 
(xi, 0) in the upper half plane y = x' > 0. Thus, integrating 

At) _ 1 

j-H{x )+H{x{t)) 
we easily obtain the expression for t(x ). This concludes the proof. □ 

Remark 3.1 The choice of ml given in the proof of Lemma 3.2 is not the 
optimal one. In fact, one could define m\ as the infimum of the ji > such 
that the level line L intersects the abscissa at some point in ]a M ,c M [. 

The next two lemmas are stated without the corresponding proofs which can 
be inferred from Lemma 3.2 via elementary considerations. 

Lemma 3.3 Let us define 

A : = sup —-^ , := sup — . 

0<s<l S 0<s<l S 

Suppose that ji e]0, mj[. Then, for every p with < po < 1 and £ with 
p < £ < 1 we have 

1 cosh-^^) < a(p ,n < J 1 . cosh-^i-). 

Note that for F like in (8) it holds that = |F'(0)|. Moreover, g - fjA > 
for fx < mj. 

Lemma 3.4 Assume (H ) and let ji> m\ (with ml as in Lemma 3.2). Then, 
for every x with 

< x < a, 

it follows that 

i— r x t ds 
lim t(x )JJI = V2 / = , 
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where Xq e ]a, 1[ snc/i £/ia£ 

P F(s) ds = 0. 



At this point, we come back to the non-autonomous system (7) and define the 
function n(t) as in [3]. Accordingly, we set 



n(t) 



rii , for < t < a, 
n , for a < t < [3, 



with 



< n < rii . 

In the sequel we treat n as a small number and rii as a large one. More 
precisely, we shall always assume < n < ml < m\ < ri\ . However, some 
further conditions on no and n\ will be required when needed. 



In order to study equations (6) and (7) we analyze the system 



(Ei) 



x =y 

y> = g X - riiF(x), 



0,1 



and follow its trajectories for i = 1 along a time interval of length a and 
for i = along a time interval of length (3 — a. Clearly, (E ) and (Ex) are 
nothing but two aspects of the previously studied system (E) with /i = n or 
H — rii , respectively. Also the energy function 8 = 8^ defined in (10) takes 
two different (but similar) forms that we denote by £ and Si , i.e., 



1 1 

£i(x, y) = £ ni (x, y) = -y 2 - g-x 2 + n^(x), i 



0,1. 



The corresponding level lines L l c (contained in the strip [0, 1] x R) are defined 
consequently. 

Let us fix a value 

ce]£ 1 (a ni ,0),0[ 
and consider the compact annular region 

Ml := {{x,y) :0<x<l,c< £ 1 (x,y) < 0} 

surrounding the point 

Pi := (a ni ,0). 
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The set Ail is invariant for the dynamical system generated by (Ei) and, 
therefore, we can use the Priifer transformation and express the solutions of 
(Ei) with initial value in Ail using polar coordinates with center at Pi . Ac- 
cordingly, we define by 9(t, z) and p(t, z) the angular and the radial coordinates 
of the solution (x(t),y(t)) of (E t ) with (x(0),y(0)) = zeM l c . 

For each point z in the set 

Af c ■= {(x,y) : a ni < x < 1, y > 0, c < £i(x,y) < 0} = M]n{ [a ni , +oo) x [0, +oo) ) 
we can fix the angular coordinate in order to have 

0(O,z)e[O,f]. 

The next result rephrases the conclusions from Lemma 3.2 and Lemma 3.4 in 
terms of the angular coordinate. Indeed, we have: 

Lemma 3.5 Assume (H Q ) and let rii > m* . Then the following properties 
hold: 

• 9(a, z) > -7r, \/zeL nAf c ; 

• for each e > 0, for each k e N and for each c G [— e, 0[ , ^/iere n* = , 
snc/i i/iai /or every ni > n*, 9(a, z) < — | — 2A;7r, \/z E L c nAf c . 

Proof. Assumption (Hq) and the choice ni > ml ensure that the level line Lq 
is homoclinic to the origin for system (Ei) (see Lemma 3.2). Thus, if we take 
z G Lq fl Af c , the solution (x(t),y(t)) departing from z moves along L in the 
clockwise sense and cannot reach the origin at any finite time. Accordingly, 
as long as we run system (Ei) we have 9(t,z) > —n, for every t > 0. This 
proves the first assertion of the lemma. On the other hand, if z G L c fl M c , 
for c < (in particular z G" Lo) we know that z belongs to a closed orbit of 
system (Ei) whose fundamental period r(z) tends to zero as p = ni grows to 
infinity (see Lemma 3.4). As a consequence, during a time interval of length 
a > 0, the trajectory will make at least L^yJ turns around Pi. From this fact, 
the second assertion of the lemma easily follows. □ 

We have now at hand almost all the needed tools to construct a domain 
containing a topological horseshoe. First, however, we need to prove another 
technical result. 

Suppose that no > is a given constant satisfying 

n < m* . 
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We consider system (Eq) and look at the orbits of such a system passing 
through a point (p , 0) with < po < 1- Recall also the definition of the 
function a in Lemma 3.2 that we call now 

0o(po,O : = / , = ds Cebo,i]- 

Jpo ^ g ^- V l)-2n^{s)-HPo)) 

The level line L° (c := £o(po,0)) crosses the vertical line x = £ exactly at 
two points and <7o(po,0 is half of the time needed to move from one of these 
intersection points to the other along the orbit of (E ) and passing through 
(Po,0). 

First of all, we observe that there is a po G ]0, a[ such that 

3- a 

<j (po, a) > — — , for every p e }0,p }. (12) 



Next, assuming (H ), we consider system (Ei) for a general n\ > ml and 
consider a level line L\ (c := ^i(po) 0)) passing through a point (p 0) 0) with 

< po < a. 

We denote by (pi,0) with 

the other intersection of such a level line with the x-axis. Note that p\ depends 
by n\ even if, for sake of simplicity in the notation, we do not make this fact 
explicit. 

We claim that there exist p e ]0, a[ and m* 2 > ml such that if we choose 

Po e]0,p o [ 

and 

n\ > m* 2 , 

then 

cro(Pi,OnJ < — ■ (13) 

To prove this claim, we need to introduce some further notation. By the sign 
conditions on F(s) and (Ho), we know that there exists a (unique) value 

b e]a,l[ 

such that 

/ F(s)ds = 0. 
Jo 
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With this position we have that 

b fM > b, V /x > ml 

and 6 M — > 6 as /x — > +00 with the function /x — > 6 M decreasing on ]m*, +00). 
Thus, if we fix a value 

/x > m* , 

we can conclude that 

b < & M < , V/x > /x. 

Hence, we can define 

77 := min -F(s) > 0. 

se[6,6 A ] 

Let us take n± — /i > p, and a point p € ]0, a[ which defines a corresponding 
point pi as described above. The proof of our claim starts now. 

As a consequence of Lemma 3.3 it will be sufficient to prove that 



Pi 



< cosh ^ - a sjg - n A^j := re, 



that is, we want to prove the inequality 

Pi>^- (14) 

K 

If ^ < a ni we are done. Hence, it is not restrictive if, from now on, we suppose 
the opposite inequality. As a consequence we have that (14) is satisfied if and 
only if the level line passing through (^, 0) has energy less than £i(po, 0). 

To this aim, we first prove that there is a lower bound So (depending on ni) 
such that 

g (b ni \2 ,b ni 



ni jrpi )>S >0. (15) 

Using the fact that £i(b ni ,0) = (by definition of b ni ), the last assertion is 
equivalent to 



We distinguish two cases, according to the fact that > b or < 6. 
In the former case, we can write 

m (^W-^t)) " f^a-^)^-^^^- |^(i-^) 

>ni?7f) -(1 - ^) := d , 
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with 5q > provided that 

„ g(K+l) 

In the latter case, we have 



= rn (^(fe ni )-^(fe)) - 1^(1 
>nmb— -(1 - ^) = 6 . 



Thus, it remains to find p such that £i(po, 0) G ] — £ , 0[ . To check this fact, 
it is equivalent to prove that 

g g 

-£i(po,0) = -n^ipo) + -pi < mrjb— (1 - 4y), 

that is, 

ni ( V + HPo) ) > f (pg + 4^ ) ■ (16) 

We define p G ]0, a[ as the solution of equation 

and take po such that 

< Po < po ■ 

For rii satisfying 

ni> - := ^V + iM) 

k(k — 1)677 

we obtain 

n 1 (,^+^(p ))>n 1?7 ^>|(a 2 + ^) 

and (16) is proved. In this manner, we have found p and our claim is proved 
for 

rii > m* 2 := max{/x*,/i}. 

We are now in position to construct a suitable domain including a topological 
horseshoe. This is discussed in the next section. 
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4 Proof of the main result 

Assume (if ) and let n(t) be like in (11), with 

< n < m* < m* 2 < n x . 

Let us fix also a number 

po e]0,p*[, with p* := mm{p ,p } 
and consider the compact annular region 

Ml := {(x,y) : <x < 1, c < S^x, y) < 0}, for c = S^po, 0). 
As in Section 3 (see, in particular, Lemma 3.5) we consider also the set 

M c ■= Ml n ( [a ni , +oo) x [0, +oo) ). 

Now we introduce a polar coordinate system with center at Pi = (a ni , 0) and 
consider the following compact subsets M' c and M" defined by 

M' c :={zeM c : e(a,z)e[- Y ,-2n]} 

and 

Ml':={zeM e : 9(a,z)e[-^,-4n]}, 

respectively. 

Notice that a point z belongs to M' c if and only if the solution (x(t) : y(t)) of 
(Pi) with (x(0),y(0)) = z e M c is such that 

(x(a), y(a)) eMln( [a ni , +oo) x (-oo, 0] ) (17) 

and x(-) has exactly two strict maximum points separated by one strict min- 
imum point along the time interval [0,a], that is the trajectory crosses the 
x— axis exactly three times. On the other hand, z belongs to M" if and only 
if the solution (x(t),y(t)) of (Pi) with (x(0),y(0)) = z e M c satisfies (17) and 
x(-) has exactly three strict maximum points separated by two strict minimum 
points along the time interval [0,a\, that is the trajectory crosses the x— axis 
exactly five times. 

Remark 4.1 For sake of simplicity in the exposition we have confined our- 
selves to the construction of only two sets M' c and M" in order to obtain at 
the end a dynamics on two symbols. We point out, however, that Lemma 3.5 
permits to define in the same manner (if one takes ri\ large enough) an arbi- 
trary number t of pairwise disjoint sets and, correspondingly, to obtain at the 
end a semiconjugation with a Bernoulli shift on £ symbols. 
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Let 7 : [0, 1] — > M 2 be a continuous mapping such that 7 crosses M c from its 
inner boundary L\ to its outer boundary L\ . Without loss of generality (if 
necessary, we can restrict 7 to a closed subinterval of its domain), we assume 
that 

£1(7(0)) = c< = £1(7(1)), 7(*) eJVc, Vte [0,1]. (18) 




Fig. 3. Example of a region Ai^ in which we have put in evidence the subset 
M c ■ The chosen parameters are the same like in Figure 2. We have taken as 
the inner boundary of M\ the level line passing through the point (0.07, 0) yield- 
ing to c ~ —0.00850436. In this case we have a ni — 0.417157 > a = 0.4 and 
£i(a ni ,Q) — —0.0936779. A path 7 crossing the inner and the outer boundaries of 
the set M c is also shown. For our argument we consider only the restriction of 7 to 
a subinterval of its domain in order to have condition (18) satisfied. 



Lemma 3.5 guarantees the existence of n* = n\ c such that for n\ > n*, 

9n 

0(a,z) < -—, V2 G L c nK, 
follows. On the other hand, 

6(a, z) > -7r, VzeL nJ\f c . 
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If we consider now the composite map 

[0,1] 3 t h-> 0(a,7(t)) G R, 

we obtain 

fl(a, 7 (0))<— 2-<-7r<fl(a, 7 (l)). 

A standard continuity argument allows to determine two subintervals [4'>4'] 
and [t[, t' 2 ] of [0, 1] with < t'[ < t' 2 ' < t[ < t' 2 < 1, such that 

7(0 G M' c , Vt G [ti, f 2 ], 7 (0 G A/" c " , Vt G 41 
and, moreover, 

0(a, 7 (4)) = - T , 0(a,7(4)) = -2vr, 

(19) 

0(a, 7 (t?)) = -y, fl(a, 7 (0) = -47r. 
Let now po be an arbitrary, but fixed number such that 

p < Po < P*- 

We denote by W the part of the plane between the level lines for energy £ of 
system (E ) passing, respectively, through (po,0) and (pi,0), where we recall 
that (pi,0) is other intersection of the level line 

£i(z,y) = £i(p o ,0), < x < 1, 

with the x— axis (see Figure 4 below). 



The intersection of M.\ with W defines two connected components 

Mlnw = AUB, 

where we have denoted by A the component contained in J\f c and by B its 
symmetric part with respect to the The sets A and B are generalized 

rectangles, that is they are homeomorphic to the unit square [0, l] 2 . 
In principle, the set A (which is the upper component of M\ fl W ) is not 
necessarily contained in Af c ■ However, such inclusion occurs when rii is suf- 
ficiently large, as we have required throughout all the paper. To check this 
claim, let us consider the intersection of an arbitrary level line 

£ (x,y) = £ (s,0), iorse[p ,Pi] 
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Fig. 4. Example of a region W C [0, 1] x R between the two level lines of energy £q 
passing, respectively, through (po , 0) and {p\ , 0) . 

with a level line 

e 1 (x,y) =e 1 (£,Q), for £g [0,p ] 

within the strip [0, 1] x R. By a symmetry argument, it is sufficient to look 
at the intersection points (x,y) with y > 0. This, in turns, yields to the 
comparison between the two functions 

(o(x, s) := ~^gs 2 + n T(s) + ^gx 2 - n T(x) 

and 

GO, : = ~2 9 ^ 2 + n ^(0 + 2 gx2 ~ n ^( x )- 

More precisely, our aim is to prove that 

&0,0 -(q(x,s) > 0, ViG [p ,o ni ], £ G [0,p ], s G [p ,Pi]- (20) 
Indeed, to check (20) it is sufficient to verify that 

CO) : = Ci(^Po) - (o(x,Po) > 0, VxG[p ,o ni ]. (21) 
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By the above positions and taking into account that T is strictly decreasing 
on [0, a], we have 

C( x ) = \a{Po 2 - Po) + n o(F(Po) - F(po)) + (no - n l )(j r (x) - F(p j) 

1 r x 

> 7,9(Po 2 ~ Po) + ( n o ~ n i) / F ( u ) du 
z Jpo 

and therefore, we can conclude that 

C(^) > \9ipQ ~ Po) > °> G [p ,Po + ], 

where po + is defined as the unique point satisfying 

rPo + 

/ F(u)du = 0, p + e]a,l[. 

Hence (21) is satisfied provided that we have 

a ni < Po + ■ 

This latter property is true for n\ sufficiently large. In fact, once that p is 
fixed, we have that p + > a is fixed too and we already observed that a M — > a 
as /i — > +00. Thus, having proved our claim, we know that the generalized 
rectangles .4. and £> look like in Fig. 5. 

We fix now an orientation for sets A and B. More precisely, given a generalized 
rectangle 71 (also called a two-dimensional cell) we select on its boundary 
two disjoint compact arcs TZf and 7Z~ , like in Section 2. Such arcs will be 
conventionally called the left and the right sides of 71. We also set 7Z~ : = 

In the case of the sets .4 and £>, we put 

^-;=^na7Wj, iT :=£n<9W. 

Accordingly, we denote by Af and A~ the two components of A~ and, simi- 
larly, by Bf and £>~ the two components of B~ (the order is not relevant). At 
last, we also define the disjoint compact sets 

Di:=^nAT c ', V 2 :=An^'. (22) 

Let ip : R 2 — > R 2 be the Poincare operator associated to system (7), that is 
z I— > C(/3, 2), where C(t) = z) is the solution of system (7) with £(0) = z. 
Observe that 

ip = ^0 ^1 , 
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Fig. 5. Example of a region M l c in which we have put in evidence the intersection 
with W consisting of two generalized rectangles A (the upper component) and B 
(the lower component). The parameters are the same like in Figure 2 and Figure 
3. Moreover, we have also chosen pQ = 0.07 and obtained po + — 0.652494. In the 
picture we have drawn the vertical lines x = a = 0.4 and x = a ni ~ 0.417157 to 
show that A C A/" c . 

where ipi is the Poincare map associated to system (Ei) for the time interval 
[0, a] and ip is the Poincare map associated to system (E ) for the time 
interval [0, (3 — a]. 

As a consequence of (19) and the fact that the annulus M.\ is invariant with 
respect to the dynamical system associated to (Ex), we have 

(V^ipt) : (A,A~) <H (B,B~), for % = 1,2. 
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On the other hand, we also have that 

(B,^ ):(B,B-)<^(A,A-). 
Thus, we can conclude that 

(V h iP) : (A,A~) <^ (A,A~), fori = 1,2. 

The results in [11] and [15], [16] guarantee the existence of a compact invariant 
set 

A C £>! U £> 2 C A 

for which is semiconjugate to a two-sided Bernoulli shift. Moreover, for 
each periodic sequence of two symbols there is a point in A which is periodic. 
In fact, all the assumptions of Lemma 2.1 are satisfied with respect to the 
Poincare map ip, the oriented rectangle (JZ.IZ^) := (A, A~) and the sets T>i 
(i = 1,2) defined in (22). 

As a final step, we want to give an interpretation of the above statement in 
terms of the solutions to equation 

x" - gx + n(t)F(x) = 0. 

Consider a two-sided sequence of two symbols (ik)kez with 

i fc e{l,2}, MkeZ. 

By the semi conjugation, we know that there is at least one point (x ,y ) = 
z G A such that 

^{z )ev lk , \fkez. 

Let x(-) (globally defined on R) be the solution of the differential equation 
with x(0) = xo and x'(0) = yo ■ For such a solution, we have that 

(x(k(3),x'(kf3))eV tk , VkeZ. (23) 

By construction, 

0<x<l, V (x, y) G T>i U T>2 ■ 
Hence, Lemma 3.1 ensures that 

< x(t) < 1, Vt G R 

and thus we have that x(-) is a solution of the original equation (6) no matter 
which kind of modification has been performed on F outside the interval [0, 1]. 

Now, we explain the meaning of (23) by analyzing the following two possibil- 
ities: 
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• {x{k(3),x'{k(3)) G V x , for some k G Z; 

• \x{k(3),x\k(3)) G V 2 , for some k G Z. 



In the former case, x(t) has precisely two strict maximum points separated by 
one strict minimum point along the time interval [(k — 1)0, (k — l)(3+a], while, 
in the latter case, x(t) has precisely three strict maximum points separated 
by two strict minimum points along the time interval [(k — l)(3, (k — Y)f3 + a]. 
In both the situations, x(-) is convex in the interval [(k — + a, k/3], with 
x'((k - 1)0 + a) < and x'(k(3) > 0. 

At last, we observe that if the sequence (ik)kez is periodic, that is ik = ik+e. 
for some £ > 1, then we can take x(-) as a £/3-periodic solution as well. 

The proof of our main theorem is complete. □ 

The proof of Theorem 1.2 relies on slight modifications of the arguments de- 
scribed above and therefore it is omitted. 

Both Theorem 1.1 and Theorem 1.2 are stable with respect to small (in the 
L 1 -norm) perturbations of the weight function n(-). 



5 Appendix 

In this section we give a proof of Lemma 3.1. Actually, we present a more 
general result which improves [20, Lemma 2.1] and may have some independent 
interest. 

We consider the second order scalar ODE 



where h : [t , t±\ x R — > K. is a Caratheodory function, that is, we assume that 
h(-,s) is measurable for all s G R, h(t,-) is continuous for almost every t G 
[to, ^i] and, for every r > there is a measurable function p r G L 1 ([to, ti], R + ) 
such that \h(t, s)\ < p r {t) for almost every t G [t , ti] and every s G [— r, r\. 
Solutions of (24) are considered in the Caratheodory sense too (cf. [7, p. 28]). 

Lemma 5.1 Suppose 



x" + h(t, x) = 



(24) 



h(t,0) = 



and, for a.e. t G [io,*i], 



h(t,s) > 0, Vs < 0, 



h(t,s) < 0, Vs > 1. 
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Let x(-) be a solution of (24) defined on [t ,ti] and such that 

< x(t ),x(ti) < 1. 

Then < x(t) < 1, Vt G [f ,*i]. Moreover, x(t) > 0, Vt G [t , *i] if h is 
(locally) lipschitzian at s = and x(t) < 1, Vt G [to,ti], ^ (locally) 
lipschitzian at s = 1. 

Proof. At first we prove that rr(£) > for all t G [to,ti]. If) by contradiction 
x(t) < for some t G [to, ti], then we can find o~o , °~\ with to < o"o < t < °"i < ti 
such that x(o" ) = x(a±) = and x{t) < for all t g]o" ,<7i[. Multiplying 
equation 

-x"(t) = h(t,x(t)) (25) 

by x(t) and then integrating on [ao,<7i], we obtain f^x'(t) 2 dt = 
J** h(t,x(t))x(t) dt < 0, a contradiction. A similar computation shows that 
v(x) < 1 for all x G [to,ti]. In fact, if by contradiction x(t) > 1 for some 
t G [to,ti], then we can find r , T\ with t < r < t < T\ < t\ such that 
x{tq) = x(ji) = 1 and x(t) > 1 for all t g]t ,Ti[. Multiplying equation 
(25) by x(t) — 1 and then integrating on [r ,ri], we obtain f^x'(t) 2 dt = 

h(t,x(t))(x(t) — 1) dt < 0, a contradiction. 
If we have minx(t) = x(t*) = 0, then also x'(t*) = since t* G]to,ti[- Hence 
x(-) is a solution of the initial value problem 

x" + /i(t, x) = 0, = x'(Q = 

and therefore, if h satisfies a Lipschitz condition at s = 0, we must have 
x(t) = 0, Vt G [t ,ti] (a contradiction with our hypotheses). 
Assume now that h satisfies a Lipschitz condition at s = 1 and maxi(t) = 
x(t*) — 1. In this (t*) = because t* G ]to, ti[ • Define also the auxiliary 

function h(t,s) := fo(t, max{l, s}) which is lipschitzian at s = 1 too and let 
y(t) be a solution of the Cauchy problem 

x" + h(t,x)=0, x(t*) = 1, x'(t*) = 

defined on a maximal interval of existence. Using the fact that h(t, s) < for 
all s G R and almost every t G [to,ti], we find that y(t) > 1 for all t ^ t* 
in the domain of y(-) and therefore y(-) is a solution of (24) as well. By the 
local uniqueness of the solutions to the Cauchy problem under consideration 
we conclude that x(-) = y(-) in a neighborhood of t* (a contradiction to 
x(t) < 1 < y(t) for t ^ t*). □ 

We remark that the same proof works for the slightly more general equation 

x" + cx' + h(t, x) = 0, 
under the same assumptions on h and for every c G R. 
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